home *** CD-ROM | disk | FTP | other *** search
/ Fritz: All Fritz / All Fritz.zip / All Fritz / FILES / PROGMISC / PCSSP.LZH / PC-SSP.ZIP / STATMISC.ZIP / POINT.FOR < prev    next >
Encoding:
Text File  |  1985-12-31  |  3.6 KB  |  133 lines
  1. C
  2. C     ..................................................................
  3. C
  4. C     SUBROUTINE POINT
  5. C
  6. C     PURPOSE
  7. C        TO COMPUTE THE POINT-BISERIAL CORRELATION COEFFICIENT
  8. C        BETWEEN TWO VARIABLES, WHEN ONE OF THE VARIABLES IS A BINARY
  9. C        VARIABLE AND ONE IS CONTINUOUS.  THIS IS A SPECIAL CASE OF
  10. C        THE PEARSON PRODUCT-MOMENT CORRELATION COEFFICIENT.
  11. C
  12. C     USAGE
  13. C        CALL POINT (N,A,B,HI,ANS,IER)
  14. C
  15. C     DESCRIPTION OF PARAMETERS
  16. C        N    - NUMBER OF OBSERVATIONS
  17. C        A    - INPUT VECTOR OF LENGTH N CONTAINING THE CONTINUOUS
  18. C          VARIABLE
  19. C        B    - INPUT VECTOR OF LENGTH N CONTAINING THE DICHOTOMOUS
  20. C          (BINARY) VARIABLE
  21. C        HI    - INPUT NUMERICAL CODE TO INDICATE THE HIGHER CATEGORY.
  22. C          ANY VALUE OF THE BINARY VARIABLE NOT LESS THAN HI WILL
  23. C          BE CLASSIFIED IN THE HIGHER OF THE TWO CATEGORIES.
  24. C        ANS - OUTPUT VECTOR OF LENGTH 9 CONTAINING THE FOLLOWING
  25. C          RESULTS
  26. C             ANS(1)- MEAN OF VARIABLE A
  27. C             ANS(2)- STANDARD DEVIATION OF VARIABLE A
  28. C             ANS(3)- NUMBER OF OBSERVATIONS IN THE HIGHER
  29. C                 CATEGORY OF VARIABLE B
  30. C             ANS(4)- NUMBER OF OBSERVATIONS IN THE LOWER
  31. C                 CATEGORY OF VARIABLE B
  32. C             ANS(5)- MEAN OF VARIABLE A FOR ONLY THOSE
  33. C                 OBSERVATIONS IN THE HIGHER CATEGORY OF
  34. C                 VARIABLE B
  35. C             ANS(6)- MEAN OF VARIABLE A FOR ONLY THOSE
  36. C                 OBSERVATIONS IN THE LOWER CATEGORY OF
  37. C                 VARIABLE B
  38. C             ANS(7)- POINT-BISERIAL CORRELATION COEFFICIENT
  39. C             ANS(8)- T-TEST FOR THE SIGNIFICANCE OF THE
  40. C                 DIFFERENCE BETWEEN THE MEANS OF VARIABLE A
  41. C                 FOR THE HIGHER AND LOWER CATEGORIES
  42. C                 RESPECTIVELY.
  43. C             ANS(9)- DEGREES OF FREEDOM FOR THE T-TEST
  44. C        IER- 1, IF ALL ELEMENTS OF B ARE NOT LESS THAN HI.
  45. C         -1, IF ALL ELEMENTS OF B ARE LESS THAN HI.
  46. C         0, OTHERWISE.    IF IER IS NON-ZERO, ANS(I), I=5,...,9,
  47. C         IS SET TO 10**75.
  48. C
  49. C     REMARKS
  50. C        THE SYMBOLS USED TO IDENTFY THE VALUES OF THE TWO CATEGORIES
  51. C        OF VARIABLE B MUST BE NUMERIC.  ALPHABETIC OR SPECIAL
  52. C        CHARACTERS CANNOT BE USED.
  53. C        THE T-TEST(ANS(8)) IS A TEST OF WHETHER THE POINT-BISERIAL
  54. C        COEFFICIENT DIFFERS SIGNIFICANTLY FROM ZERO.
  55. C
  56. C     SUBROUTINES AND FUNCTION SUBPROGRAMS REQUIRED
  57. C        NONE
  58. C
  59. C     METHOD
  60. C        REFER TO P. HORST, 'PSYCHOLOGICAL MEASUREMENT AND
  61. C        PREDICTION', P. 91 (WADSWORTH, 1966).
  62. C
  63. C     ..................................................................
  64. C
  65.       SUBROUTINE POINT (N,A,B,HI,ANS,IER)
  66. C
  67.       DIMENSION A(1),B(1),ANS(1)
  68. C
  69. C     COMPUTE MEAN AND STANDARD DEVIATION OF VARIABLE A
  70. C
  71.       IER=0
  72.       SUM=0.0
  73.       SUM2=0.0
  74.       DO 10 I=1,N
  75.       SUM=SUM+A(I)
  76.    10 SUM2=SUM2+A(I)*A(I)
  77.       FN=N
  78.       ANS(1)=SUM/FN
  79.       ANS(2)=(SUM2-ANS(1)*SUM)/(FN-1.0)
  80.       ANS(2)= SQRT(ANS(2))
  81. C
  82. C     FIND NUMBERS OF CASES IN THE HIGHER AND LOWER CATEGORIES
  83. C
  84.       P=0.0
  85.       SUM=0.0
  86.       SUM2=0.0
  87.       DO 30 I=1,N
  88.       IF(B(I)-HI) 20, 25, 25
  89.    20 SUM2=SUM2+A(I)
  90.       GO TO 30
  91.    25 P=P+1.0
  92.       SUM=SUM+A(I)
  93.    30 CONTINUE
  94. C
  95.       Q=FN-P
  96.       ANS(3)=P
  97.       ANS(4)=Q
  98.       IF (P) 35,35,40
  99.    35 IER=-1
  100.       GO TO 50
  101.    40 ANS(5)=SUM/P
  102.       IF (Q) 45,45,60
  103.    45 IER=1
  104.    50 DO 55 I=5,9
  105.    55 ANS(I)=1.E38
  106.       GO TO 65
  107.    60 ANS(6)=SUM2/Q
  108. C
  109. C     COMPUTE THE POINT-BISERIAL CORRELATION
  110. C
  111.       R=((ANS(5)-ANS(1))/ANS(2))* SQRT(P/Q)
  112.       ANS(7)=R
  113. C
  114. C     COMPUTE T RATIO USED TO TEST THE HYPOTHESIS OF ZERO CORRELATION
  115. C
  116.       T=R* SQRT((FN-2.0)/(1.0-R*R))
  117.       ANS(8)=T
  118. C
  119. C     COMPUTE DEGREES OF FREEDOM
  120. C
  121.       ANS(9)=FN-2
  122. C
  123.    65 RETURN
  124.       END
  125. (I)
  126.    30 CONTINUE
  127. C
  128.       Q=FN-P
  129.       ANS(3)=P
  130.       ANS(4)=Q
  131.       IF (P) 35,35,40
  132.    35 IER=-1
  133.       GO